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Abstract 

We discuss the effective potential of the conformal factor in the effective average 
action approach to Quantum Einstein Gravity. Without invoking any truncation or 
other approximations we show that if the theory has has a non-Gaussian ultraviolet 
fixed point and is asymptotically safe the potential has a characteristic behavior near 
the origin. We argue that this characteristic behavior has already been observed in 
numerical simulations within the Causal Dynamical Triangulation approach. 



1 Introduction 



One of the problems one faces when one tries to construct a fundamental theory 
of quantum gravity is the absence of any experimental data or observations that neces- 
sarily would have to be explained by the theory [1-4]. In this situation it is particularly 
important to compare the predictions of a priori different candidate theories and to find 
out whether they are perhaps just different formulations of the same underlying theory 
or whether they really belong to different "universality classes". In this way one can at 
least narrow down the set of independent possibilities among which the experiment must 
decide in the end. Guided by the experience with Yang-Mills theory we would expect 
that in particular the comparison of continuum and lattice approaches should be very 
instructive and fruitful. On the side of the continuum approaches, recently a lot of efforts 
went into the exploration of the asymptotic safety scenario [5-28] in the formulation based 
upon the gravitational average action. It aims at defining a microscopic quantum field 
theory of gravity in terms of a complete, i.e., infinitely extended renormalization group 
(RG) trajectory on the theory space of diffeomorphism invariant functionals of the metric. 
The limit of an infinite ultraviolet (UV) cutoff is taken by arranging this trajectory to 
approach a non-Gaussian fixed point (NGFP) at large scales (k — > oo). This NGFP of the 
effective average action is not only instrumental in constructing the quantum field theory 
by dictating how all generalized couplings must "run" when the UV regulator scale is sent 
to infinity, it also determines the physical properties of the resulting regulator-free theory 
at large physical scales, the behavior of propagators at large momenta, for instance. We 
refer to this quantum field theory of the metric, defined in the continuum by means of 
the effective average action, as Quantum Einstein Gravity (QEG). 

In Yang-Mills theory, the asymptotic scaling behavior that controls the approach 
to the continuum of the corresponding lattice gauge theory can be read off from the 
corresponding effective average action [29-32]. It is therefore plausible to ask whether 
there is an analogous relationship between QEG and gravity on a (dynamical) lattice. 
A first comparison (of critical indices) has been attempted in [9] where the results from 
the average action were confronted with Monte Carlo data obtained in Regge calculus 
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simulations [33]. The results were found to be consistent, albeit within considerable error 
bars. 

A different statistical-mechanics based framework for quantum gravity that has been 
developed during the past years is the Causal Dynamical Triangulation (CDT) approach 
[47-50]. It provides a background-independent, nonperturbative definition of quantum 
gravity in which the causal (i.e., Lorentzian) structure of spacetime plays a crucial role. 
Only causally well-behaved geometries contribute to the path integral. It is regularized by 
summing over a particular class of triangulated, piecewise flat geometries. All contributing 
spacetimes are foliated by a kind of proper-time and can be obtained by gluing together 4- 
simplices in a way that respects this foliation. The slices of constant proper-time consist 
of spacelike tetrahedra forming a 3-dimensional piecewise flat manifold. Typically its 
topology is assumed to be that of a three-sphere S 3 . Each configuration contributing 
to the regularized path integral can be rotated to Euclidean signature, and the rotated 
partition function can be computed by means of standard Monte Carlo methods. It was 
found [50] that the CDT approach can generate almost classical 4-dimensional universes 
on which only small quantum fluctuations are superimposed. 

In ref. [48] an effective action for the Robertson- Walker scale factor a(t) of these 
universes has been extracted from the Monte Carlo data. In the present paper we are 
going to compare this effective action to the corresponding prediction of QEG. We shall 
demonstrate that the latter reflects the non-Gaussian fixed point, the crucial prerequisite 
of asymptotic safety, by a characteristic "finger print" at small scale factors. The very 
same finger print has actually been observed already in the CDT simulations. 

We interprete these results as an indication pointing in the direction that QEG 
and the CDT approach could possibly describe the same underlying theory in different 
languages, i.e. that they are two representatives of the same "universality class". 

In fact, in ref. [8,23] another point of contact between CDT and QEG had been 
found already: They agree on the microscopic spectral dimension of macroscopically 4- 
dimensional spacetimes; in either case one finds the somewhat surprising result ds = 2. 

The remaining sections of this paper are organized as follows. In Section [2] we 
analyze the effective potential of the conformal factor in QEG and discuss in particular 
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the impact the NGFP has on it. In Section [3] we make contact with the CDT simulations, 
and in Section @] we present the conclusions. 

2 The Effective Potential of the Conformal Factor 
in QEG 

In this section we discuss the standard effective potential (i. e., the one with vanish- 
ing infrared cutoff, k — 0) for the conformal factor of metrics on maximally symmetric 
spacetimes with the topology of a d- dimensional sphere S d . The starting point is the exact 
gravitational effective average action [6] along some RG trajectory, Tk[g^v, g^u], and the 
related reduced functional Tk[g^] = ^k[g^u, g^v\- (The ghost arguments are set to zero 
and are not indicated explicitly.) The latter functional is assumed to have a representation 
of the form 

ffcW = ^2u a (k) i a [g^\ (2.1) 

a 

where {/ a [(yV^]} is an infinite set of local and nonlocal "basis" functionals, invariant un- 
der diffeomorphisms acting on g^ u , and the w Q 's are the corresponding running coupling 
constants. We denote their canonical mass dimensions by d a = [u a ]. Hence, since l\. is 
dimensionless, [I a ] = —d a . The dimensionless running couplings are defined by 

u a (k) = k~ da u a (k) 

so that we may rewrite (12. ip as 

ffcfo*,] = k da Ia[9»u] (2.2) 

a 

Up to now the metric argument g^ v was completely general. At this point we 
specialize for metrics on S d , with a variable radius 0. We parametrize them as 

9 iw = 4> 2 9txv (2.3) 
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where g^ is the metric on the round S d with unit radius, and the conformal factor <fi is 
position independent. Hence is a metric on a round sphere with radius <fi. We shall 
denote the volume of the unit-S d by a d = J d d x y/g = 2^ d+1 ^ 2 /T((d + l)/2). 

We use conventions such that the coordinates x M are dimensionless and <fi has the 
dimension of a length. Hence [# M „] = —2, and g^ u is dimensionless, [g^u] = 0. 

Without having made any approximation so far, the effective average potential for 
the conformal factor, £4(0), by definition, obtains by inserting the special argument (12.31) 
into 

U k (<j>) J d d x Jg~ = T k [ 9ilv = 4> 2 g, u ] (2.4) 
In terms of the expansion (12.21) we have the exact representation 

M4>) = cr^J^Mk) k^Utfg^] (2.5) 

a 

or, more explicitly, 

UM = a^J^Mk) (k<P) da Ia[9„u] (2.6) 

a 

To obtain equation (12.61) we exploited that I a [(p 2 g^u] = 4> da Ia[9[i.A which holds true since 
I a has dimension — d a . (This relation can be regarded the definition of the canonical mass 
dimension.) 

Eq. ( I2.6P makes it manifest that if we know a complete RG trajectory {u a (k),0 < 
k < oo} we can deduce the exact running potential from it, and in particular its A; — ^ 
limit, the standard effective potential [7 e fr(0) = £4=o(0)- Usually we are not in the com- 
fortable situation of knowing trajectories exactly; nevertheless certain important proper- 
ties of U e s(4>) can be deduced on general grounds. For this purpose we shall employ the 
following decoupling argument which is standard in the average action context [29,31]. 

The basic observation is that the true, i.e. dimensionful coupling constants u a {k) 
have a significant running with k only as long as the number of field modes integrated 
out actually depends on k. If there are competing physical cutoff scales such as masses or 
field amplitudes the running with k stops once k becomes smaller than the physical cutoff 
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scales. (See Appendix C.3 of [28] for an example.) In the case at hand this situation 
is realized in a particularly transparent way. The quantum metric is expanded in terms 
of eigenfunctions of the covariant (tensor) Laplacian D 2 of the metric g^. This metric 
corresponds to a sphere of radius 0; hence all eigenvalues of the Laplacian are discrete 
multiples of 1/0 2 . As a result, when k has become as small as k ~ 1/0, the bulk of 
eigenvalues is integrated out, and the a a 's no longer change much when k is lowered even 
further. Therefore we can approximate 

UM) = f4 =o (0) « £/*=i/*($ (2.7) 

In order to make the approximation (12. 7p strictly valid we have to be slightly more 
specific about the precise definition of Uk(<f))- The above argument could be spoiled by 
zero modes of D 2 . Therefore we define and Uk in terms of a functional integral over the 
fluctuation modes of the metric with a non-zero eigenvalue of D 2 only. As a result, the 
actual partition function would obtain by a final integration over the zero modes which 
is not performed here. The only zero mode relevant in the case at hand is the fluctuation 
of the conformal factor. When we compare our results to the CDT simulations later on 
it will be important to keep in mind that U e s(4>) has the interpretation of an effective 
potential in which the conformal fluctuation has not yet been integrated out since it is 
essentially this quantity that is provided by the Monte Carlo simulations. 

Eq. (12.71) has a simple intuitive interpretation in terms of coarse graining: By low- 
ering k below 1/0 one tries to "average" field configurations over a volume that would be 
larger than the volume of the whole universe. As this is not possible, the running stops. 
Note that the S topology enters here; the finite volume of the sphere is crucial. 

With the approximation (12.71) we obtain the following two equivalent representations 
of U e g(4>) in terms of the dimensionless and dimensionful running couplings, respectively: 

um) = ^^Mr 1 )/^] (2.8) 

a 

C4ff(0) = a^^air 1 )^^) (2.9) 

a 

As an application of these representations we consider two special cases. 
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Let us assume the RG trajectory under consideration has a classical regime between 
the scales k\ and k%, meaning that u a (k) ~ const = u^ lass for k\ < k < ki- Then (12.91) 
implies that for k^ 1 < <t> < , approximately, 

UM) = afY,^ 1 ^]^ (2.10) 

a 

As expected, this potential has a nontrivial ^-dependence governed by the classical cou- 
plings M^ lass . 

Next let us explore the consequences which a non- Gaussian fixed point has for the 
effective potential. We assume that the dimensionless couplings u a (k) approach fixed 
point values w* for k — ■> oo. More precisely, we make the approximation u a (k) « u* a 
for k > M with M the lower boundary of the asymptotic scaling regime. Then the 
representation (12.81) tells us that £7 e ff(0) = crj" 1 J2 a u a Ia[gnA if ~ M -1 . Obviously this 
potential is completely independent of (p: 

U eS ((f)) = const for all < M _1 (2.11) 

In typical applications (see below), M equals the Planck mass mpi = £p\ so that U c s is 
constant for <fi < l-p\. 

Eq. ( 12. lip is one of our main results. It shows that the existence of an ultraviolet 
fixed point has a characteristic impact on the effective potential of the conformal factor: 
Regardless of all details of the RG trajectory, the potential is completely flat for small 
<f). The interpretation of this result is that for <fi < M _1 the cost of energy (Euclidean 
action) of a sphere with radius (f) does not depend on <p. Spheres of any radius smaller 
than M~ l are on an equal footing. This is exactly the kind of fractal-like behavior and 
scale invariance one would expect near the NGFP [8,10]. 

We emphasize that except for the decoupling relation (12.71) no approximation went 
into the derivation of this result. It is an exact consequence of the assumed asymptotic 
safety, the existence of a NGFP governing the short distance behavior. Neither has the 
theory space been truncated nor have any fields been excluded from the quantization 
(such as in conformally reduced gravity the transverse tensors, for instance, cf. [21,22]). 
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On the basis of the above general argument we cannot predict how precisely, or how 
quickly the effective potential flattens when we approach the origin. However, we expect 
that its derivative dll e g/d(j) vanishes at = 0. This has an important physical implication. 
In general, possible vacuum states of the system (the "universe" ) can be found from the 
effective field equation 5T k=0 / Sg^ = 0. More specifically, S^-type groundstate candidates 
have a radius 0o at which (dU e g / d(ft) (0o) = 0. Thus we see that thanks to the NGFP a 
vanishing radius 0o = has become a vacuum candidate, the derivative of U e g vanishes 
there. (To qualify as the true vacuum it should be the global minimum.) Hence the 
universe has an at least metastable stationary state with = 0, i.e. a state with a 
vanishing metric expectation value (g^ u ) = 0. In this state gravity is in a phase of 
unbroken diffeomorphism invariance, see ref. [22] for a discussion of this phase in the 
context of asymptotic safety. 

Let us finally illustrate the above discussion in the familiar setting of the Einstein- 
Hilbert truncation [6] in d — 4 which is defined by the ansatz 

f k W = ~Y^G~ k j d ' x Vg{R{g) - 2A fc ) (2.12) 

Inserting (12.31) with = 0(x) we obtain 

r ^/-] = jJuTa / d 4 x - \r u d^d v 4> - 2 + (2-13) 



For x-independent only the potential term survives, with 

3 / ,9,9 1 



( - k 2 <p 2 + ^X(k) fc 4 4 



Ang(k) V 6 

If A(A;) > 0, the case we shall always consider in the following, Uk{4>) has a minimum at 
a nonzero radius given by 

0o(AO = VWW) (2-15) 
This is exactly the radius of the S 4 which solves the ordinary Einstein equation following 



from the action ( l2.12ll R In the second line of (12. 14ft we employed the dimensionless 



because this space is maximally symmetric, by Palais' theorem [34], inserting the ansatz <? M „ = (j^g^i, 
commutes with deriving the critical point. 
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Newton constant g(k) = k 2 G(k) and cosmological constant X(k) = A(k)/k 2 . So there are 
the following two equivalent ways of writing the effective potential: 



u eS (<t>) = I 

3 



1 .2,1 A(0 



+ 



(2.16) 



The RG trajectories of the Einstein-Hilbert truncation have been investigated and 
classified in [9]. Here we can concentrate on those with a positive cosmological constant, 
those of "Type Ilia", because this class of trajectories seems to be the one which is behind 
the Monte Carlo results we are going to discuss. Important regimes along a Type Ilia 
trajectory include 

The NGFP regime: g(k) « g*, X(k) « A* for k > M. 

The k 4 regime: G(k) « const, A(&) oc k 4 for k-r < k < M, where fcx is the "turning 
point" scale at which f3\ vanishes. 

The classical regime: G(k) ~ const = G, A(k) ps const = A for fcterm ^ k < k^ where 
fcterm is the scale at which the Einstein-Hilbert truncation breaks down and the trajectory 
terminates at a singularity^. 

If one defines the classical Planck mass and length by mp\ = £p^ = G^ 1 ^ 2 one finds that, 
approximately, M « mpi. (For further details see [9,40,44]; see in particular Fig. 4 
of [44].) 

In the /c 4 -regime, when k decreases, the cosmological constant quickly becomes 
smaller proportional to k 4 , and the radius of the sphere "on shell", <po{k), increases pro- 
portional to 1/k 2 . 

If the underlying RG trajectory of QEG is of Type Ilia then U e s(4>) is constant 
in the NGFP regime < Ipi, and it equals the classical potential for A;^ 1 < < fc t " e J m . 
Note that our ignorance about the infrared end of the trajectory entails that we have no 
information about the effective potential for large values of (p. The intermediate /c 4 -regime 



2 If one tentatively matches the trajectory against the observed values of G and A one finds that 
fcr ~ 10 _30 TOpi, corresponding to k^ 1 w 10 _3 cm, and fcterm ~ 10 -60 mpi ss Hq so that k^~ eim equals about 
the present Hubble radius [44], [40]. 
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of the trajectory gives rise to a behavior 



U eS {(p) oc (-0 2 + const) for k^ 1 < < £ Ph 



(2.18) 



In the above discussion we tacitly assumed that the trajectory is such that M 2 « 
mpj 3> A; otherwise no classical regime would exist. 



Figure 1: The effective potential for the Type Ilia trajectory discussed in the text. The 
dashed line represents the potential t/ c i ass with the same values of G and A, but all 
quantum effects neglected. 

A qualitative sketch of the resulting U e g is shown in Fig. [TJ It is compared there to 
the classical potential L^gs which would obtain if G and A had no ^-dependence at all. 
The crucial difference between the two is the almost constant U e g at small <fi. This regime 
is a pure quantum gravity effect, directly related to the existence of a NGFP. Quantum 
mechanically, but not classically, the universe can be stationary at small values of 0, at 
least at (f> = 0. 

As a consequence of our assumption A m| b the U e s = const regime ends at a 
radius (f> ~ £p\ which is smaller than the classical "on-shell" radius 4>o = \/3/A. The 
actual "size of the universe" corresponds to a scale in the classical regime of the RG 
trajectory therefore. 



t/eff 
A 




9 



In the region where the quantum effects modify [/class most strongly the term oc cf 2 
is the dominant one. We can therefore say that the key effect behind the flattening of 
the potential near the origin is the running of Newton's constant. Its consequence for the 
shape of U e g can be understood as the result of the "RG improvement" [35-46] 

k - wh^f = 7, (2 ' 19) 

with G(k) = g*/k 2 , as appropriate near the NGFP. 



3 Making Contact with CDT Simulations 

The CDT approach [47-50] defines a discrete version of the Wick rotated quantum- 
gravitational proper-time propagator 

Gl G [g 3 (0),g 3 (t)] = J Vg E e~ s ^ (3.1) 

Here Se is the Euclidean Einstein-Hilbert action, and the integration is over all 4-dimensio- 
nal Euclidean geometries g E of topology S 3 x [0, 1], each with proper-time running from 
to t, and with prescribed spatial boundary geometries #3(0) and ^(t), respectively. In the 
simulations reported in [48], for technical reasons, periodic rather than fixed boundary 
conditions have been used so that the topology of the spacetimes summed over is S 3 x S 1 
rather than S 3 x [0,1]. (Furthermore, the simulations were done at constant 4- volume 
V4 rather than constant A; the corresponding propagator is related to (13.11) by a Laplace 
t r ansf or mat ion . ) 

It is instructive to visualize the typical, statistically representative 4-geometries 
contributing to the path integral. They are characterized by a function V^s), < s < t, 
where V^s) is the 3- volume of the spatial S 3 at proper-time s. If t is large enough, 
a "typical universe" has long epochs with a very small V3 at early and late times (the 
"stalk") and in between a region with a large Vs(s), see Fig. 1 of ref. [48]. 

It has been shown [48] that the dynamics of these "universes" is well reproduced by 
a minisuperspace effective action for rotated Robertson- Walker metrics 

ds 2 = dt 2 + a 2 (t) dQl (3.2) 
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where df2g is the line element of the unit 3-sphere so that Vs(s) oc a 3 (s). It reads 

For large a the potential V is 

V(a) = -a + ^Aa 3 = V cl {a) (3.4) 

The action (13.31) with (13.41) is, up to an overall minus sign, what one obtains when one 
inserts (13. 2p into the Einstein-Hilbert action. (In simulations with fixed V4 the constant 
A is a Lagrange multiplier to be fixed such that J" * ds V^s) = V4.) For small a, the 
(— a)-term in H3.4[) is to be replaced by a function of V3 whose derivative at a = behaves 
as V^, v > 0. A convenient parametrization, valid for all a, is 




Figure 2: The effective potential V(a) of eq. (13. 5p obtained from the Monte Carlo simula- 
tions (solid line). The dashed line shows the potential V c \(a) from the classical Einstein- 
Hilbert action. (From Loll et al. [48]) 

This potential is sketched in Fig. [2] and compared to V c \(a). While the latter has a 
negative slope at the origin, V^(0) = — 1, the potential (13.51) with the corrected small-a 
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behavior has vanishing derivative there: V'(0) = 0. This quantum induced modification 
is important since (13.51) allows for a classically stable solution a(t) = which explains the 
stalk observed in the computer simulations [48]. 

To summarize: According to the Monte Carlo simulations, the effective potential 
of the Robertson- Walker scale factor agrees with the classical one for large a, but differs 
from V c \ at small a in a characteristic way: V(a) becomes essentially constant (zero) for 
a — > 0. In the parametrization (13.51) the classical V(a) ~ —a ist turned into V(a) ~ — |a 3 
by the quantum effects, which amounts to an almost flat a-dependence for a<l. 

Comparing Figs. [Hand [2] makes it obvious that in both cases the dominant quantum 
gravity effects have an analogous impact on the effective potentials U c g(<f)) and V(a), 
respectively: Their slope at the origin , while strictly negative classically, vanishes thanks 
to the quantum effects, and for small 0, or a, the potential is essentially flat, rendering a 
state with a = or = stationary. 

It is plausible to assume that we are dealing with the same phenomenon in both 
cases. In fact, upon introducing the conformal time rj{t) = J*dt' /a(t') the line element 
(13. 2p assumes a form analogous to (12.31) . 

ds 2 = <j)(r ] ) 2 [dr ] 2 + dn 2 3 ] (3.6) 

with the conformal factor <j)(r)) = a(t(rj)). Since <fi and a differ only by a time reparametriza- 
tion, which is irrelevant here, the potentials U e g((f)) and V(a) are almost the same object. 
In particular we defined U e g(4>) in terms of a functional integral (or the corresponding flow 
equation) which does not include the conformal zero mode, i. e. fluctuations which merely 
change the radius of the S 4 . Likewise its CDT counterpart V(a) results from integrating 
out all modes other than the spatially constant global scale. 

Thus we see that within QEG the non-Gaussian fixed point predicts a vanishing slope 
of Ues at (ft = 0, and upon identifying U e s with V(a) this feature is actually observed in 
the CDT simulations: They all generate universes with a "stalk", and according to the 
analysis in [48] this stalk owes its existence to the vanishing derivative of V(a) at a = 0. 
In a certain sense, the stalk seems to realize the "exotic" phase of gravity when the 
metric has zero expectation value and the group of diffeomorphisms is not broken by the 
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groundstate (see also section 6 of [22]). 

Another minor difference between the QEG and CDT setting, respectively, is that 
g^y is a metric on S* 4 , while drf + &VL\ refers to S 3 x [0, 1] or S 3 x S 1 . However, we do 
not expect such global issues to cause qualitative changes for small conformal factors^. 
(It is also interesting that with the stalk removed the computer-generated universes [50] 
are almost perfect 4-spheres S A .) 

4 Summary 

We discussed the effective potential of the conformal factor both from a QEG and a 
CDT perspective. We demonstrated that if QEG is asymptotically safe then it gives rise 
to a potential which becomes flat for cf) — > 0, allowing for a phase of gravity with vanishing 
metric expectation value. The argument assumes the existence of an underlying UV fixed 
point, but is exact otherwise. We argued that the potential V(a) "measured" in Monte 
Carlo simulations within the CDT approach does indeed reflect the predicted behavior 
at small scale (or conformal) factors. It is therefore plausible to assume that already the 
presently available CDT simulations, indirectly, have seen evidence for the same NGFP 
that has been found in QEG by truncated flow equations. This hints at the possibility 
that QEG and CDT might be related at a deeper level. 



Acknowledgement: M.R. would like to thank J. Ambj0rn, H. Hamber, R. Loll, and 
R. Williams for helpful discussions. 



3 The motivation for taking g^ v to be a S 4 metric in Section [2] is that only in the case of maximal 
symmetry it is easy to identify the decoupling scale in terms of geometric data. 
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